Intraplanar Magnetic Excitations in NaiCo0 2 : 
An Inelastic Neutron Study 
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Inelastic neutron scattering measurements mapping the in-plane magnetic interactions of 
NaiCoCb reveal dispersive excitations at points above an energy gap E p = 11.5(5) meV at the 
superstructural Bragg reflections. The excitations are highly damped, broadening with increasing 
energy, and disappear at hoo ^35 meV, a strong indication that the magnetism is itinerant. Tilting 
into the ac plane reduces the value of E^ by 25%, suggesting that the dispersion along c is significant 
and the magnetic correlations are three-dimensional, as seen at the higher doping levels. 

PACS numbers: 61.12.-q, 61. 12. Ex, 67.57.Lm, 71.27.+a, 75.30.Fv, 72.10.Di 



The discovery of superconductivity in hydrated sodium 
cobaltate in 2003 1 has generated a large volume of study 
in a relatively short space of time. The anhydrous par- 
ent phase is interesting in its own right; the system, at 
x>0.65, is a Curie- Weiss metal with magnetic ordering 
below 22 K and the signature of magnetism, assigned 
to a spin density wave (SDW), depends sensivitely on 
the level of sodium doping via its influence on the Co 
valence state. Recent neutron investigations of the SDW 
phase 2,3 demonstrated that the magnetic ordering can be 
ascribed to an A-type antiferromagnetic (A-AFM) struc- 
ture. In this system, however, only 1 — x of the Co sites 
carry spin, so that the simple Hamiltonian used to de- 
scribe the excitation modes would only be valid if phase 
segregation into magnetic and nonmagnetic regions is as- 
sumed. This raises the question of whether the system 
really can be characterised by a model with localised spin 
ordering or whether an itinerant model with small charge 
disproportionation is more appropriate. A new regime 
emerges at x = |, with different magnetic ordering. At 
Tm = 88 K, long-ranged static AFM ordering occurs 
which has been clearly identified in neutron diffraction 
experiments 4 . At Tmi = 53 K, the material undergoes 
a metal to insulator transition (MIT). Further to this, 
there exists long-ranged ordering of the Na cations 5 into 
quasi-onedimensional chains along a, enabling the crystal 
structure to be reclassified into an orthorhombic super- 
cell of dimensions (2a x \/^a x c). In this crystal setting, 
there are now two Co sites. This led to speculation that 
the above transitions were the results of charge ordering 
(CO) in the C0O2 planes, probably due to the influence 
of charge modulations in the Na layer. More recently, 
however, the notion of a charge ordered Mott-like tran- 
sition has been disputed. NMR measurements 6 did not 
deduce any large difference in charge state between the 
two sites. Crystal structure determination by powder 
neutron diffraction 7 also failed to find a significant dif- 
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ference between the Co-0 hybrization on the two sites. 
A clear picture of the role of spin and orbital ordering 
in NaiCo02 has yet to be defined. Here we present the 
first measurements of inelastic neutron scattering on this 
system to probe the nature of the interactions in the ab 
planes. A dispersive mode, magnetic in origin, is ob- 
served; the high damping and large velocity of the excita- 
tion indicates the possibility that the AFM ground state 
is itinerant. Preliminary investigations indicate that the 
system is also highly dispersive along c. 

The samples were prepared by the floating zone 
method as described previously. To confirm the stoi- 
chiometry of each sample, given that the magnetic tran- 
sitions at Tmi and Tm occur only at x = ~, magnetic 
susceptibility (x) measurements were carried out in a 
Quantum Design MPMS-5S superconducting quantum 
interference device (SQUID) magnetometer on a range of 
samples; those selected for the experiments below had x 
similar to previously reported measurements^**^. The re- 
sults indicated the sample was free from magnetic impu- 
rities (CoO or C03O4). The heat capacity (C) measure- 
ments were carried out by a two-tau relaxation method 
in a Quantum Design PPMS. Inelastic neutron scattering 
measurements were taken on the MAPS spectrometer at 
ISIS (RAL) and the triple axis spectrometer (TAS) 2T1 
at the Laboratoire Leon Brillouin (LLB) (CEA-CNRS). 
A mosaic of six crystals, each of mass ~0.5g, was used for 
the experiment at ISIS. The samples were loaded into a 
CCR and cooled to 10 K. An incident neutron energy of 
100 meV was selected and long counting times (typically 
48 hours) were adopted to build up sufficient statistics in 
(Q,E) space over the range of interest. Just one of these 
crystals was selected for the triple axis measurement in 
order to reduce the overall mosaic spread and optimise 
resolution. The crystal was cooled to 2.8 K in an orange 
cryostat with the c axis aligned vertically so that the Co 
moments were in the scattering plane. 

Magnetic susceptibility and heat capacity measure- 
ments are shown in figure ^ Both the magnetic ordering 
and the metal-insulator transitions are clearly visible as 
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FIG. 1: Magnetic susceptibility (R a b = 5 T)and heat capacity 
of Nai C0O2. The two transitions at Tm and Tmj are clearly 
visible in the magnetisation. On the contrary, the heat capac- 
ity data is featureless. The magnetic structure, first proposed 
by Choy et al. 9 is presented in the inset, both the hexagonal 
(shaded blue) and orthorhombic (pink) unit cells are shown. 

shoulders in the magnetisation. It is interesting to note 
that x does not exhibit the Curie- Weiss-like character 
as seen in the SDW material: there is no signature of 
charge locality from the macroscopic measurements. The 
magnetic structure as measured by Gasparovic et al£ is 
presented in the inset. The symmetry breaking imposed 
by the orthorhombic supercell produces two separate Co 
sites, one of which is magnetically ordered (Co 4+ S=|) 
and the other nonmagnetic (Co 3+ S=0). In terms of the 
new unit cell, the system can be thought of as a simple G- 
type antiferromagnet with a magnetic propagation vector 
of (l|T). The dynamics of the magnetism may be more 
complicated than suggested by this localised moment or- 
dering picture; fluctuations on the nonordered rows of Co 
may affect the spin wave dispersion. The fact that the 
MIT occurs 35 K below Tm would suggest the transport 
carriers reside mainly on the second Co site. The opening 
of a gap on the part of the Fermi surface (FS) related to 
band containing these carriers would readily explain the 
onset of the insulating transition. 

In order to confirm the nature of the magnetic order- 
ing proposed by Gasparovic et aim the crystals were first 
rotated by 30° to the incoming neutron beam in order to 
mimic the experimental setup as used by Boothroyd et 
al™, where in-plane ferromagnetic (FM) ordering for an 
x = 0.7 crystal was visible as scattering around the ori- 
gin, corresponding to a zero 2D magnetic ordering vector. 
No such scattering exists at this sodium level confirm- 
ing the absence of FM coupling. The crystal was then 
rotated back to have k^ || c. Inelastic scattering inten- 
sity is clearly visible (figure EJ) at points above (^00) in 
the parent hexagonal lattice. All six structural Bragg 
peaks are also visible in the elastic channel, which is 
at odds with a true long-ranged orthorhomic superstruc- 
ture, where there should only be two. However, the inte- 
grated intensity of pairs of peaks diametrically opposite 
each other match, suggesting that the crystals exhibit 
twinning (rather tripling) in the basal plane and are not 



just aligned along different equivalent direction choices of 
(100). S(Q,fkj) maps were produced by slicing the data 
across different energy transfer levels. The small size of 
the sample resulted in a low signal to background ratio; 
the background itself is a product of multiple and instru- 
mental scattering and increases with |Q|. The energy 
resolution is dependent on the opening time of chopper 
used to monochromate the beam; this was relaxed to 
collect good statistics but ensured that elastic scattering 
extended up to 10 meV. As a result the background ob- 
scures the magnetic scattering at |Q| >0.75. Magnetic 
excitations should be visible as a cone of scattered in- 
tensity above the magnetic zone centre. There are no 
in-plane magnetic Bragg reflections 4 and so, to further 
investigate the origin of the scattering seen here, mea- 
surements were completed at the TAS 2T1 at LLB. 

At a temperature of 2.8 K, a clear peak is identifi- 
able in both the inelastic and elastic channels (a non- 
magnetic superstructural peak resulting from Na order- 
ing 5 ) at Q=(^00). In order to access high enough energy 
transfers to examine the dispersion, it was necessary to 
move to the next Brillioun zone. The magnitude of the 
change in intensity between the two sites is comparable 
to that expected from the decrease in the magnetic form 
factor of Co 4+ with the increase in |Q|. A further indica- 
tion that the signal is a result of 3D magnetic ordering is 
that the peak is destroyed by increasing the temperature 
to ~2Tm (figure |3J). Similar sets of Q-scans were made 
at increasing energy transfers and are shown in figure 21 
The excitation is gapped, with no scattered intensity visi- 
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FIG. 2: Inelastic scattering in the (hkO) plane summed over 
energy transfers of 10-20 meV. Scattering is visible at points 
corresponding to [(100)] (|00) in the [orthorhombic] hexagonal 
[super-] cell, although twinning reproduces the hexagonal sym- 
metry. An overlay of the hexagonal reciprocal lattice is shown. 
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FIG. 3: Constant energy scans through (ffO) at fuo — 15 
meV. The superstructural Bragg peak at this Q, and there- 
fore the Na ordering, persists to high temperatures, whilst 
magnetic Bragg peaks at (j>j>lodd) disappear at Tm 4 . The 
inelastic signal has disappeared at twice the ordering temper- 
ature Tm , with an increased background corresponding to an 
enhancement of the incoherent and multi-phonon scattering. 



ble below ~11.5 meV. The peaks are also overly damped, 
increasing in width with energy, whilst the overall dis- 
persion is very steep. The data can be fitted to a linear 



dispersion huj = y A 2 + (DQ) where A is the value of 
the energy gap and D the spin wave stiffness constant, 
more commonly referred to as the spin wave velocity in 
AFM systems. A fit to the data gives A = 11.5(5) meV 
and D = 510(20) meV A and is plotted as the solid line 
in panel B of figure El 

To examine the excitation gap further, scans at con- 
stant Q were made with increasing energy transfer. The 
gap was observed at both (|-00) and (4-|0) (figure [£)) 
at E^= 11.5 meV. The onset of scattering is sharp, and 
decays with increasing fwj. As seen at the higher dop- 
ing levels, there exists an optical phonon (OP) at ~20 
meV which dominates the signal in the Q-scans at that 
Huj. Since OP modes are relatively flat, it was possible 
to pick up the signal from the magnetic excitation once 
more at higher hco. In order to correct for the additional 
phonon component to the scattering, another energy scan 
was taken at an offset position from the excitation at 
(0.594 0.4 0), the exact values of (hkO) being chosen to 
set the analyzer at the same 20 position as the previ- 
ous scan. The difference curve then removes the incoher- 
ent [phonon] scattering below [above] the gap. Finally, the 
sample was tilted slightly into the (100)-(001) scattering 
plane and the same energy scan repeated for Q=(j>0j>) 
(not shown). The value of is shifted to ~9 meV sug- 
gesting that the dispersion along c is significant, a feature 
in accord with the 3D magnetic correlations reported for 
x >0.7— . Q = 0) is not a magnetic zone centre, no 
magnetic scattering exists at this point. In order to de- 
termine whether the excitation includes a spin wave gap, 
information of the dispersion along c is necessary. 



These data exhibit features that are common to in- 
tinerant magnetic systems such as Cr— or MnsSiA^, in 
which the highly damped dispersion, at the position of 
a magnetic Bragg peak, is seen to vary linearly in Q 
up to energies of ~2 kgT/v and then rise vertically in 
a "chimney" of scattering. The dispersion here, how- 
ever, remains linear in Q until the signal disappears at 
~4 keT/v- The most obvious method for determining the 
proportion of localised to itinerant spins in a FM material 
is to calculate the ratio of the high temperature param- 
agnetic moment and the low temperature saturated mo- 
ment, equal to one in completely localised magnetic sys- 
tems and greater than one for itinerant materials. This 
method is not applicable here; \ is n °t Curie- Weiss-like 
and the material has AFM ordering. For a localised spin 
system, the energy of the excitation at the zone boundary 
should be of the order of k^T^. An extroplation of the 
data collected in this study predicts these two energies to 
be different by an order of magnitude. If the system ex- 
hibits low dimensionality, the energy scale characterising 
the excitation will relate to the strength of interactions 
along the spin chains, thus underestimating k^Tjvf. How- 
ever the loss of the excitation into a broadened continuum 
of decay modes, as witnessed here, indicates strongly that 
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FIG. 4: A: Constant energy transfer scans at Q=(f j>0) taken 
along (100) at 2.8 K. The datasets have been offset by -150, - 
200 and -250 counts respectively for increasing energy transfer 
cuts. B: The measured spin wave dispersion along (hjjO) (the 
variation of the spin wave damping F with energy is shown in 
the inset); the excitation appears at ~11.5 meV. 



4 



600 



(Vi Vi 0) 

<y 2 +s Vi-h oy 

net counts 




15 20 25 

Tico (meV) 

FIG. 5: Energy scans at (j-§0). A background was taken 
at a position of (-| + J ~ — JO) where J ^0.1 was chosen so 
that the spectrometer was aligned with the same 2$ positions 
as the on-peak measurements. The subtracted data (black 
open squares) is plotted; the line corresponds to the convolu- 
tion of the experimental resolution ellipsoid and the intensity 
expected for an ideal gapped system (E^ = 11.5 meV). 



the material is itinerant in character. 

In the light of these results, the previously proposed 
model of a charge ordered ground state should be recon- 
sidered. No evidence of either the Tm or Tmi transition 
can be seen in the specific heat. Recent NMR studies 6 
have suggested that both transitions could be attributed 
to spin density waves i.e. FS nesting effects, in which 
case Tmi would involve only modest modifications to 
the spin order, producing little measureable magnetic 
entropy. Alternatively, as mentioned above, the mag- 
netism in NaiCo02 may exhibit low dimensionality; in- 
plane correlations build slowly with decreasing tempera- 
ture before locking in-phase along c at Tm- The phase di- 
agram of sodium cobaltate is characterised by the effects 



of strong electron correlations. At x>0.5, some of the Co 
states may be localised by these interactions causing the 
Curie- Weiss-like magnetic behaviour and spin wave ex- 
citations, whilst the material retains metallic conductiv- 
ity. Below half doping the electron correlation effects are 
weaker and the system is Pauli paramagnetic. Therefore 
there must exist a cross-over between these two regimes, 
a special state realised at x=0.5 in which the magnetism 
is largely descibed by the itinerant electron picture but 
electron correlation effects play a role in the formation 
of the MIT. It should be noted that the insulating tran- 
sition is anisotropic; the resisitivity at low temperatures 
along c is 25 (]cm and just 15 mllcm along afrAA, a value 
comparable to that of the "metallic" SDW system 8 , and 
whilst NaiCo02 is considered a low temperature insu- 
lator in dc, optical conductivity measurements-^ reveal 
the MIT is only visible in the low frequency data. The 
mechanism that drives the MIT is clearly not yet under- 
stood and further investigation into the electronic band 
structure at this doping level is needed. 

In summary, the intraplanar magnetic excitations 
whose zone centres coexist with structural Bragg reflec- 
tions have been studied. The excitation has a gap of 
11.5(5) meV. The localized spin model does not pro- 
vide physically realistic values for the magnon velocity 
or the energy transfer at the zone boundary and is not 
consistent with the broadening and diminishing scatter- 
ing more reminiscent to the transformation between spin 
wave and Stoner excitations as seen in the itinerant mag- 
nets. The excitation modes measured in this study in- 
dicate that the system is at least partially itinerant and 
that a model containing zero spin Co 3+ and spin-half 
Co 4+ is not appropriate. 
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